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ABSTRACT 


Approximate formulas for the inband radar cross section of arrays with 
parallel feeds are presented. To obtain the formulas. multiple reflections are 
neglected, and devices of the same type are assumed to have identical electrical 
performance. 

The approximate results were compared to the results obtained using a 
scattering matrix formulation. Both methods were in agreement in predicting RCS 
lobe positions, levels, and behavior with scanning. The advantages of the 
approximate method are its computational efficiency and its flexibility in handling 


an arbitrary number of coupler levels. 
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I. INTRODUCTION 


The design of radar-stealthy platforms has become an important engineering 
problem. The principles of radar stealth have been well known for several decades, 
but only recently have technological advances allowed practical implementation of 
these principles. Some examples of stealthy platforms are the SR-71 spyplane. the 
F-117A fighter. and the B-2 bomber. All of them have low radar cross section 
(RCS), low infrared (IR) emissions by control of heat sources, and low microwave 
emissions by using "quiet" radar and communications. 

The RCS of future military platforms will be lowered significantly by the use 
of shaping and materials selection. Consequently, attention is now being focused 
on the onboard sensors in order to ensure that their signatures do not become 
predominant, Of particular interest are wideband phased arrays, which can possibly 
operate at the same frequencies as those of an illuminating threat radar. Thus a 
high performance array must not only meet the system antenna operating 
requirements (gain. sidelobe level. etc.) but also the RCS requirements. It is 
essential that any technique applied to minimize the radar cross section does not 
seriously degrade the primary operation of the antenna system. 

Solid state microelectronic devices now permit integrated transmit/receive or 


receive-only elements to be collocated with array apertures. Self-calibrating and 


adaptive “smart skin" conformal arrays are now practical. Increased efficiency. 
reliability, and reduced cost are the benefits of this technology. This technology 
makes phased arrays more appealing than other antenna types. 

In characterizing and, subsequently, minimizing the RCS of phased arrays, 
reliable computation of scattering is crucial. It is necessary to consider both inband 
and out-of-band threat frequencies because the scattering characteristics of an array 
are distinctly different in the two frequency regions. If the incident wave is in the 
array's operating band, the radiating elements are well matched and threat signal 
can penetrate into the feed and be reflected at internal mismatches and junctions. 
Even for a well-matched array, there can be a large number of scattering sources 
that add constructively under some conditions. This effect depends on the type of 
feed and the devices incorporated therein. In some cases. these reflections can 
significantly modify the RCS. 

In this thesis, the inband RCS of a phased array with a parallel (corporate) 
feed. as shown in Figure 1, is examined. The main objectives are: 


1. to find approximate equations for the inband RCS of phased arrays with 
parallel feed networks, 


to 


to compare the approximate solution with a rigorous solution based on 
scattering parameters, and 


Lvs) 


to determine the RCS behavior for various feed parameters for both linear 
and two-dimensional arrays. 
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| Array Antenna with Parallel Feeds 


Figure 1. Typica 


In general. antenna scattering analysis is very difficult. Hansen [Ref. 1] 
describes two antenna scattering modes. The first one is the antenna or radiation 
mode. which is determined by the radiation properties of the antenna. and vanishes 
when the antenna is conjugate matched to its radiation impedance. The second is 
the structural mode. which is generated from currents induced on the antenna 
surfaces. The two modes are not easily identifiable. in particular, when the array 
is installed on a platform. In this thesis only the antenna mode is examined. which 
is the dominant RCS component for a phased array in its operating band. 

Chapter I provides the theoretical background for array antenna RCS analysis 
and includes the derivation of the approximate formulas. It presents information 
about the correspondence between RCS lobes and the location of periodic scattering 
sources within the feed network. Chapter III describes the formation of the so- 
called rigorous solution. This method includes multiple reflections between devices 
in the feed whereas the approximate solution only considers the first reflection and 
neglects higher order reflections. Chapter IV deals with the comparison and 
analysis of results. Finally, Chapter V concludes with a discussion on the benefits. 
concerns. and recommendations that can be applied to the analysis of low 


probability of intercept radar and communications systems. 


lf. THEORETICAL BACKGROUND 


Antennas with identical amplitude and phase patterns can differ in the way 
they scatter. The possible importance of this point in connection with evaluation 
of antennas designed to have identical patterns was first examined by R. H. Dicke 
{Ref. 2], who proposed "to see what can be done in the way of differentiating 


between a good and bad antenna on the basis of scattering." 


A. DEFINITION OF RCS 

The assigning of a radar cross section to a given target, whether it is an 
aircraft or a ship, is based on the fact that this object functions as an antenna, The 
back-scattering cross section may therefore be interpreted as a measure of the 
antenna currents excited on such an object. Radar cross section is a measure of 
power scattered in a given direction when a target is illuminated by an incident 
wave. (An older term for RCS is echo area.) Mathematically, the RCS is defined 


as 


E (6)/? 
o(8,) = lim4nR? E@: (1) 


Row |E (8) |? 


where R is the distance from the target to the observation point (receiver), E,(8) is 
the scattered electric field in the direction of the receiver. and £,(0,) is the incident 
electric field (assumed to be a plane wave), The term monostatic means that the 
transmitter and receiver are co-located with respect to the target, @ = 8,. RCS has 
units of square meters. 

In general, RCS is a function of the angular orientation and shape of the 
scattering target. as well as frequency and polarization of the transmitter and 


receiver. Some typical values of RCS are listed in Table 1. 


Table |. RCS of Common Targets 





Insects 


Fighter Aircraft 


Bomber Aircraft 




















The scattering characteristics for all targets fall into three natural regions. 
The first one is the low-frequency or Ravleigh region. defined by &L<<1 (L is the 
characteristic length of the target and k = 27/A is the wave number). The second 
region is the Mie or resonance region. where AL is on the order of unity. The third 
is the high-frequency or optical region. where AL>>!. An array of length L can 
fall into either of the three categories depending on the threat radar's wavelength 


relative to the antenna's operating wavelength. 


B. SCATTERING FUNDAMENTALS 

Antenna scattering has been the subject of study since 1950, however there 
is very little published on the subject in the open literature. Until the mid 1980's, 
work has concentrated on the analysis of low gain antennas. Recently. high gain 
antennas have received attention because it is anticipated that they will be in 
widespread use on future low RCS platforms. 

For an antenna subjected to an incident electromagnetic wave. the scattered 
field has traditionally been divided into two modes: the antenna mode (£,) and 


the structural mode (£,) [Ref. 1]. The total scattered field (£) is the sum 


iZ, ] 
J” hh-E) £—|r (2) 
AR, R 





EZ‘) = E(Zj) +E, = EZ) + 





where 
Z,= R, + jX, is the radiation impedance 
fi = effective height of the antenna 
Zo = 377Q 
E, = incident field 
The structural mode arises from currents which are induced on the antenna 
and the surrounding structure when the terminating load is equal to the complex 
conjugate of the antenna impedance. The antenna mode results when the induced 
current delivered to the antenna feed point is reflected and then reradiated. The 
antenna mode is proportional to the gain of the antenna for a given direction and 


a modified reflection coefficient, which is equal to 


cS (3) 





where Z, is the impedance across the antenna terminals. 

When the antenna is conjugate matched. the modified reflection coefficient 
(r) goes to zero and the antenna mode vanishes. Usually, this is achieved when the 
antenna impedance has a real value and the terminals are connected to a matched 
transmitter and/or receiver. Consequently, there is no reflection at the junction 


between the antenna and the transmission line. However. Z, for an element in an 


array is dependent on the angle of arrival of the incident wave because of mutual 
coupling. It is practically impossible to force + to zero at every angle 
simultaneously. 

If there is a mismatch as defined by the modified reflection coefficient. then 
antenna mode reflections occur. The amount of antenna mode scattering depends 
on the gain of the antenna in a particular direction. Thus, if the antenna does not 
have high gain outside of its operating band. then it will not have significant 
antenna mode scattering even though there is a large modified reflection 
coefficient. 

Kahn and Kurss [Ref. 3] have demonstrated that for a large class of antennas 
terminated by matched receivers, the scattered power is generally greater than the 
absorbed power. equality being attained for minimum-scattering antennas. This 
result has frequently been interpreted to mean that no conjugate-matched antenna 
can absorb more than it scatters. This implies that gain must be lowered to reduce 
RCS. However, Green [Ref. 4] has shown that an antenna can absorb more than 
it scatters if the gain in the back direction exceeds its gain in the forward direction. 


This is the case for all high performance phased array antennas. 


C. SCATTERING CHARACTERISTICS OF PHASED ARRAY 
A single dipole provides low directivity. To increase antenna size. and hence 


directivity. a collection of elements can be arranged and interconnected to form an 


array. The basic element can be an aperture or slot, horn, microstrip patch, spiral, 
or dipole depending on the application. In order to provide very directive patterns, 
the fields from the array elements must interfere constructively (add) in the desired 
directions and interfere destructively (cancel each other) in the remaining space. 
The RCS of an array antenna can be decomposed into the components 
described in the previous section: the antenna mode and the structural mode. The 
relative importance of the two terms will depend primarily on the threat frequency. 
As shown in Figure 2 the frequency domain is separated into five bands. For a 
well-designed antenna, the RCS in the operating band should be low because most 
of the incident energy is delivered to the antenna load. However, even though the 
individual reflections from the antenna are small, a large array can have have tens 
of thousands of such sources. Thus the RCS can achieve significant levels under 


some conditions. 


Low out-of- Transition In-band Transition 
band band band 


Figure 2. Antenna Frequency Band 





High out-of- 
band 





In this thesis, only threat signals in the operating band of the antenna will be 
considered. In this case. the wave penetrates into the feed network and is reflected 


at internal junctions and devices. The total RCS is determined by the vector sum 


of the individual scattered fields that return to the aperture and reradiate. These 
scattering sources are depicted in Figure | for a parallel feed network. They 


include: 


1. aperture. 7, 


lo 


phase shifter inputs. 7, 


inputs of the first level of couplers, r, 


Goa 


4. loads at the sum and difference arms of the first level of couplers. To te 


ta 


sum and difference arm loads for higher levels of couplers. etc. 


Typical sources of scattering are the following: 


1. Devices that are not perfectly matched due to physical limitations 
(manufacturing and assembly errors, variation in the electrical properties of 
material, etc.) 


to 


Surface roughness (assembly tolerances. discontinuities, distortions. etc.) 
3. Errors built into the antenna (e.g.. quantization errors from subarraying) 


4. Edge effects. 


The inband antenna mode RCS is obtained from equations (1) and (2) and 


summing over all array elements 


4nA; 


9 
“ 





N i 
7,008 2 |r, (9) 2 m 


net 


o(8,~) = 


11 


where 
T,(8.¢) = total reflected signal returned to the aperture for element n 
when the wave is incident from the (0,@) direction 
k= k(ut + vy + w2) = wave vector 
u = sinBcos@ 
v = sinOsing 
w = cos0 
d, = position vector to element 7 
. = effective area of an element = /°Z,/4R, 


FE 


nora 


= normalized element scattering pattern 


To arrive at (4) identical elements have been assumed; the variation in mutual 
coupling near the array edges has been neglected. This allows the RCS to be 
separated into an array factor and an element factor, just as in the radiation case. 

Evaluating (4) requires the total reflected field at each element. A rigorous 
solution must employ a network matrix formulation such as scattering parameters. 
If multiple reflections within the feed can be neglected, an approximate solution 


can be obtained by tracing signals through the feed and back to the aperture. 


D. 


APPROXIMATE METHOD 


In this section approximate RCS formulas are derived. The following 


assumptions are made: 


lo 


All the devices of the same type are assumed to have identical electrical 
characteristics. That means that all the radiating elements have the same 
reflection coefficient r, and the same transmission coefficient ¢. None of 
the elements is ideal, because each reflection coefficient r, is not equal to 
zero. By the same token. all phase shifters have a reflection coefficient r,, 
etc. 


All couplers are represented by magic tees. which implies equal power 
splitting. (This is not a low sidelobe feed). 


In the operating frequency band, all feed devices are well matched and 
therefore higher order reflections are neglected (r << 1). 


Only scattering from the aperture, phase shifter inputs, coupler inputs, and 
the sum and difference arms of the first and second levels of couplers are 
considered. Couplers in higher levels of the network are assumed to be 
perfectly matched. 


Lossless devices are assumed for simplicity, which implies 


In|? + |e\? = 1 (5) 


for a device where r is the reflection coefficient. and ¢ is the transmission 
coefficient. 


Identical aperture elements with a Lambertian scattering pattern (cos*6). 


Edge effects are not included. 


8. 


Assuming that only one scattering source dominates at any given angle, the 
coherent sum of the scattered signals is represented by a noncoherent sum 


[byt Bj te # EY eB ee, ann 1B? (6) 


where E, is the reflected signal of the n® element. Thus, the total RCS can 


be expressed as 6 = 6, + 6, + Oy + Oy + Oy + Og t 


Random errors are neglected since they only contribute to an average RCS 
level. 


Parallel feeds are suited to rectangular element arrangements, and therefore 


linear and rectangular array geometries will be studied. Array quantities are 


defined in Figure 3. Note that: 


th 


wa 


For the case of linear arrays, all elements are aligned along the x-axis and 
equally spaced, d. The z-axis is broadside to the array. 


For the case of two dimensional planar arrays all elements are in the xy- 
plane, uniformly spaced with dimensions d, and d,, and numbers of elements 


N, and N,. 


There is only a 8 polarized incident field. (For linearly polarized elements 
in the xy plane. this gives rise to the cos@ scattering pattern.) 


The phase shift per element introduced by the phase shifter in Figure | is 
A. Furthermore, the phase shifters are reciprocal. 
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Figure 3. Two-Dimensional Array Geometry 


The reflection sources for an eight element linear phased array with a parallel 
feed are shown in Figure 1. The incident plane wave at an angle @ arrives at each 
radiating element, which is the first scattering source encountered with reflection 
coefficient r,. The transmitted signal, which is determined by the transmission 
coefficient ¢,, proceeds to the phase shifter. Again, if the phase shifter is not 
matched to the transmission line, a reflected signal returns to the aperture, 


depending on the reflection coefficient r,. A portion of this reflected field is 


reflected again by the radiating element since it is assumed to be reciprocal. This 
is a second-order reflection. and will be neglected in the calculation of RCS. The 
portion of the signal not reflected at the phase shifter is transmitted to the first level 
of couplers, and so on. 

Assuming that each phase shifter has an insertion phase, the signal that 
propagates through it will encounter a phase shift depending on the antenna beam 
scan angle @,. For a linear phase progression. the transmission coefficient for the 


phase shifter at element # is 


ee ie! (7) 

where x, = Adsin®, for the linear array. 

Figure 4 shows the first magic tee in the array with the fourth port loaded. 
(It can represent any three-port power divider as well.) Usually mismatches exist 
in both side arms, and in the sum and difference arms. The mismatches result in 
reflections back to the aperture, with reflection coefficients of r, at ports 2 and 3. 
rs at port 1, and r, at port 4. The angle A is the signal phase at the coupler port 
3 relative to the signal phase at the coupler port 2. and includes all of the insertion 
phases of the devices between the couplers of the first level and the aperture. It 
also includes any space path delay relative to the origin for the incident wave. If 
the phases of all reflection coefficients (except the phase shifter) are zero. then 


A = @ - % where a = kdsind. 
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Figure 4. Reflections at the Ports of a Magic Tee. 


The magic tee shown combines signals from the first two radiating elements. 
Similarly, the next magic tee combines signals from the third and fourth elements 
(24 at port 3 and 3A at port 2). The portion of the signal that enters the sum 
arm propagates down to the second level of couplers where some reflection occurs. 


with the remaining signal transmitted to the third level of couplers. and so on. For 


a parallel feed network, the number of radiating elements is always of the form 2” 
where m is the number of levels of couplers. 

The signals reflected from the magic tee can be determined from its scattering 
matrix. First assume that the magic tee is perfectly matched and therefore has the 


following scattering matrix: 


01 0 
echo o 1 (8) 

2/100 -1 

Id 1-1 0 


Referring to Figure 4. the input signals at ports 2 and 3 are £, and E, 
respectively. Both signals have unit amplitude but differ in phase by A radians. 
The subscripts on E refer to the element numbers to which tee's side arms are 


connected. In the sum arm the combined signal is 


ee ee 

! 2 2 
A (9) 

IZ 
= j2e ° ese). 
2 
Similarly, for the difference arm the combined signal is 
Ei SE, By:= v2 _ y2 ia 


4 


aS 
= ape? sin). 


The next magic tee in the array combines the input signals E, and E, from 
elements 3 and 4. Again. both signals have unit amplitude but differ in phase by 
A (because £, has a phase of 2A and £, has a phase of 3A) 

E. = E, +E, = v2 ,ja + v2 ,u3a 
= 3 4 ,) 


(1) 
A 


and 


= Jaa 
E, = E, -E, = ~j/Ze? sin(S) . (12) 
In general. the signals in the sum and difference arms are 
ga 
Ey = y2e : cos efln-2a (13a) 
: 5) 
and 
if A 
E, = -jy2e sin 5 ef (13b) 


for n=1],2...N/2. 
The reflected signals trom the sum and difference arms of the first level of 


couplers have the form 


Ey = Ey ty (14a) 


and 


r.. (14b) 


Therefore the total reflected signals returning at the side arm inputs are 

















.A 
, , ' LF 
E, = Ey 8, + Ey Sy, = reco 5 | jrgsi 5 Crs, (15) 
el A AV 2 
E, = Ey 83, + Ey Sy, = reo + irssin( $})e os (16) 
A A)| 3 
E; = racoy{ >) ~jrgsiol 5) e > eA | (17) 
A A\] (3 
Eye2 racos{ >) + jrysio( 3) Bele ™ (18) 





As expected, these equations indicate that the RCS depends on the relative phases 
of the signals entering the side arms. 
The total scattered signal due to first level of couplers is obtained by 


summing all reflections returned to the aperture 


20 


ES = E+E, + E,+E,+..+ Ey. (19) 
The phase shifter delay and path delay in the direction of the observer must be 


included. Summing the terms in (19) vields [Ref. 5 ] 


oe 2( A\/ sin(NA) sin?( 4 \{_sin(VA) | 
EY = Ty COS 4) - 5 20 
5 sin(2A) fa 7 caer sina) a 





This can be interpreted as the sum of returns from the sum arm plus returns from 


the difference arms. Applying the definition of RCS gives 








o. = SRA? 22 coef AY[_sin(va) F 
E, ae ened || ee (21a) 
— sin(2A) 
2 
and 
‘ 4nA? 2. (>) sin(NA) F 
a > A a (21b 
: N sinc2a)| 


The portion of the signal not reflected by the first level of couplers is 
transmitted down to the second level where the entire process is repeated. The 


total signal returned to the aperture due to reflections at the second level of 


couplers is 


Ee = 4/Bcos'| 5 ‘[rycos*A - r,sin2a}| SiN) sige (22) 
2 sin(4A) 


Again, this is a sum of RCS contributions from the sum and difference arms, 





2 ‘ 2 
0, = 4nA rs cos'(3} ‘costA ° _Sin(NA) 





3 a2 (23a) 
— sin(4A) 
4 
_ 4nA2 2 fA). 4, .[ sin(NA) P 
Fa, = rr rx COS (3) sin"A Wi a (23b) 
poe) 


The coupler RCS contribution contains three factors. The first one is 
(4n47r°/X’), the RCS of a reflector of area A reduced by the reflection coefficient 
, The squared terms in large brackets are an array factor for the couplers. 
Finally, the remaining factors are equivalent to an element factor for the coupler 
sum or difference arms. 

For linearly polarized radiating elements along the x-axis and a 0-polarized 


incident wave, the element factor is 


F__ (6,9) = cos0. (24) 


norm 


This factor can be lumped with 4 to form a projected area. The physical area is 


ined 
i) 


related to the number of elements and the spacing. Thus for a linear array 


A = Ndlcos@ (25) 


where 

N is the number of elements 

d is the spacing between elements 

/ is the effective height of the element. 
For every RCS contribution. the first and the last factor of (4) remain the same. but 
the second factor is different. These are designated as array factors even though 
they include the subarray factors for the couplers. The array factor for the 
radiation element contribution is: 

AF, = ran): (26) 
Nsin(a) 

The array factor for the phase shifter contribution is: 


AF =r tla (27) 
Pp PT | Ngin(a) 


The arrav factor for the side arms of the first level of couplers becomes: 


2,2. ( sin(NA 
AF, = t,t; rae) (28) 


The array factor for the sum arms of the first level of couplers becomes: 


sin(NA) 


N 


ae of A 
AF, = t,t ry COs () 
ee) 


(29) 








/ 


which occurs in (21a). The array factor for the difference arms of the first level 


of couplers becomes: 


AF, = ¢tr s'( 5) _sin(VA) | 


rtyry 30 
: > sin(2A) oy 


which also occurs in (21b). 


The array factor for the sum arms of the second level of couplers becomes: 


AF, = iets cost 5) cos*(A) { SN(NA) 
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which is related to (23a). The array factor for the difference arms of the second 


level of couplers becomes: 


AF, = ce rs cos!) sin?(A) ass) 
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which is also related to (23b). By combining (4) and (6). the total monostatic 
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inband antenna mode RCS is 


4nA? 


2 2 2 2 2 2 2 
a AF, + AF, + AF, + AFs, + AF, + AF, +AFs). (3) 


Oo = 





Consider a of two-dimensional planar array in the xy-plane with a rectangular 
grid as shown in Figure 3. The array has N, by N, elements spaced d, by d,. 
Parallel feeds are used to combine signals for elements along the x-axis. Thus the 
results for the linear array can be applied directly to the two-dimensional array by 
choosing A properly and multiplying all terms by a y-direction array factor. If 
0, and @, are the scanned antenna beam angles. then let 


Xo, = kd, sin®, cos , (34a) 


and 


Xoy = kd, sin®, sine, (34b) 


be the interelement phases to scan the antenna beam. Now A in the linear array 


formulas is replaced by ¢ where 
C, = kd,u - Xp, - (35) 


All scattering terms for the linear array must be multiplied by an array factor for 


the y-dimension. For scattering sources ahead of the phase shifters 
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ap. ~ (Sn, kd,)) (36) 
> (N,sinkd,») 


and for scattering sources behind the phase shifters 


_ (80,6) (37) 
»  (Nysin(Z,) }’ 
where 
Gy. = REY > Koy (38) 


Thus. for a two-dimensional array equations (26) through (32) become 


ee sin(N,kd,u) sin(N, kd, v) (39) 
"| N,sin(kd,u) N,sin(kd,v) 


AF =r? 
P oe. Famed 





sin(N kd,u) sin(N kd, vy) (40) 
sited 





AF, = 


c rit, 


2.2_ (sin(N,¢,)\(sin(N,¢,) (41) 
N,sin(¢,) } | N,sin({,) } ” 


c,|{ sin, ¢,) | fsindy, 6) 
AF, = iit coe © , ae faa , (42) 
ri sin(2¢,)| 7 7” 
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a al S)f sin, 6) \/ sin, Z) 
APs 7 Ta sa §) ered a) (43) 
—sin(2¢,) vee ¥ 
5 
AF, = ‘ t rs co cos*(¢,) sits. $) ae ; (44) 
sina ¢,) i 





APs = bly cost sin?(¢,)| ne s*o? ee i) 


N\sin(¢, (45) 

—sin(4f)}\ > (oy) 
4 

Finally. the total RCS is 


7 4nA? 
a2 





(AF, + AF, + AF) + AFs + AF; +AFz,+AF,) (46) 
where 


A = N,N,d,d,cos® . (47) 


Ill. RCS ANALYSIS FOR THE RIGOROUS SOLUTION 


In the approximate method, the higher order reflections are neglected because 
they vary as 7°, 7°, etc., where r << 1. A rigorous solution based on a scattering 
matrix formulation contains the effects of multiple reflections. Furthermore, it is 
possible to combine scattering matrices with the method of moments to include 
interactions between the feed and aperture as well as reflections inside the feed. 
This method solves the problem rigorously by obtaining an antenna impedance 
matrix that describes the electrical characteristics of the antenna surfaces and feed 
network. The method of moments (MM) impedance matrix is combined with the 
feed scattering matrix and continuity equations that relate the MM expansion 
coefficients to the feed signals as described in [Ref. 6]. 

In this thesis a variation of the rigorous method in [Ref. 6] is used, where the 
radiating element is represented by a simple two-port device with reflection 
coefficient r,. This eliminates the method of moments portion of the antenna 
matrix in the rigorous solution of [Ref. 6], and the problem reduces to a pure 
scattering matrix one. This allows direct comparison of the rigorous solution with 


the approximate solution, because the radiating elements are modelled the same. 
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A typical scattering matrix equivalent network with eight elements is shown 
in Figure 5. If only two levels of couplers are modelled. then only a four-element 
array need be considered. Reflections trom the third level of couplers are zero and 
therefore adjacent four-element subarrays are decoupled. The four-element array 
problem can be solved by the scattering matrix formulation in [Ref. 6]. 

In the rigorous method. the number of radiating elements is assumed to be a 


multiple of 4. The incident plane wave field incident on element / : 


So etre. PS has (48) 


t 


The N elements of the array can be divided into N/4 subarrays as shown in Figure 


6, where the a coefficients are used to obtain the scattered field F° 


N 
en jkdu(n -1) 
E ye 
rel 
NA 4 ™ 
7 j2kdu(n - 14 jkdu(i-1)4 
ye J ay eltout- | 
ne i=] 


Equation (49) has been obtained using the fact that the four-element subarrays are 


decoupled 


Gy = Ug = Ay = -- (50) 


The total RCS is obtained from (1). 
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Figure 5. Typical Linear Array with N Elements Divided into N/4 Subarrays. 
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4 - element linear array 
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Figure 6. Scattering Matrix Equivalent Network when the Radiating 
Element Match is Constant with Angle 

For each four-element subarray, there are eight 2-port devices (four radiating 
elements plus four phase shifters) that yield 16 equations, and three 4-port devices 
(two magic tees in the first level plus one magic tee in the second level) that yield 
12 equations. There is a total of 28 equations and 28 unknowns (28 a's) for the 
four 4-element subarray. 

The important characteristics of the rigorous method are the following: 


1. Interactions between all devices in the feed and radiating elements are 
included; 


to 


The transmission line phase between devices is included. Thus the effect 
of changing line lengths can be modelled. 
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IV. RCS OF LINEAR AND TWO-DIMENSIONAL ARRAYS 


In order to simplify the analysis and reduce the demand for computer time. 
only linear arrays have been examined using the rigorous method. Thus a 
comparison of the two methods is only given for linear arrays. RCS contour plots 
for two-dimensional arrays were also obtained to illustrate RCS behavior with 


beam scanning. The MATLAB programs are included in the Appendices. 


A. RCS DATA FOR THE APPROXIMATE SOLUTION 


1. Linear array 

For the approximate method, RCS pattern data for linear arrays was 
computed. The linear array program computes the RCS per squared wavelength 
in dB for any monostatic angle 8, number of radiating elements N, and scanned 
angle 8. The contributions from each scattering source can be broken out 
individually if desired. 

The inband RCS of linear arrays with 16, 64, and 128 elements for 
Q, = 0 degrees (no scanning) is shown in Figures 7 through 9. All reflection 


coefficients for feed devices and the radiating elements are equal to 0.2 
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Figure 7. Inband RCS of a Linear Array with a Parallel Feed for N=16, 
6, = 0° (Approximate Method). 
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In Band RCS of a linear array with parallel feed N= 64 and thetas=0 
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Figure 8. Inband RCS of a Linear Array with a Parallel Feed for N=64, 
8,= 0° (Approximate Method). 
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In Band RCS of a linear array with parallel feed N= 128 and thetas=0 
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Figure 9. Inband RCS of a Linear Array with a Parallel Feed for N=128, 
6, = 0° (Approximate Method). 
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For all linear array calculations the spacing is d = 4/2 and the effective height (or 
length) of the elements in the y direction is /= A/2. 

Referring to the figures, note that the number of major lobes (spikes) is 
the same in all cases. For a small number of elements (e.g., N = 16) the lobes are 
not as well defined because they are broader and lower. The RCS contributions 
from first and second levels of couplers are shown in Figures 10 and 11, 
respectively. The lobe spacing in the RCS pattern for the first level of couplers is 
determined by the physical spacing of the couplers (2d). This can be generalized 
for higher levels also. For instance, the effective spacing of the second level of 
couplers is 4d. Therefore, as more levels of couplers are added to the feed, more 
lobes appear between already existing lobes. 

Figure 12 illustrates the effect of beam scanning. Assuming that the phase 
shifters are reciprocal devices, the lobes associated with mismatches behind the 
phase shifters scan with the antenna beam because of the factor xy). The large 
lobe. at 9 = 45 degrees in Figure 12 is due to in-phase addition of the scattered 
signals passing through the phase shifters, as expected. Note that the specular lobe 
at 8 = 0 degrees does not scan. The high lobes near + 85° are due to Bragg 


diffraction (the RCS equivalent of grating lobes). 
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Figure 10. Contribution from First Level of Couplers. 
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In Band RCS of a linear array with Parallel feed N= 64 and thetas=0 
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Figure 11. Contribution from Second Level of Couplers. 
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Figure 12. Inband RCS of a Linear Scanning Array with a Parallel Feed for 
N = 64, 6, = 45° (Approximate Method). 


39 


2. Two dimensional array 

A second MATLAB program, which is shown in Appendix B, was used 
to generate contour plots of inband RCS for two-dimensional arrays. For 
simplicity, only square arrays are examined ( N, = N, and d, = d, = A/2) for a 
specified scan angle (6,, 9,). Again it is assumed that all reflection coefficients 
are 0.2. RCS contours are shown in Figures 13 through 16. They are plotted in 
direction cosine space (uw and v) and the contours enclose spatial regions of CS 
above a specified level. This level is chosen to be 10 dB when N, = N,,= 16. and 
20 dB when N, = N, = 64 elements. 

The only difference between the arrays of Figures 13 and 14 is the 
number of elements: 16 by 16 versus 64 by 64. Both plots are symmetric about 
the horizontal axis at v= 0. High RCS is present along the principal planes of the 
array and then drops off away from these axes because of the separable product in 
the array factors. It can be seen that the lobes are narrower and higher as the 
number of elements increase. The specular lobe appears at the center of each plot. 

Figures 15 and 16 show the RCS of scanned arrays. The first one has 
8, = 45 degrees and ~, = 0 degrees, and the second 0, = 45 degrees and ~, = 45 
degrees. Both plots are symmetric around the axis v= 0. As in the linear array 
case, the lobes originating from reflections behind the phase shifters scan along 


with the antenna beam. 
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In Band RCS of a two-dimensional array - Contour plot 
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Figure 13. Inband RCS of a Two-Dimensional Array with a Parallel Feed for 
N, = N, = 16. 8, = 0°, o, = 0°. Contour Plot at 10 dB Level 
(Approximate Method). 
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In Band RCS of a two-dimensional array - Contour plot 





Figure 14. Inband RCS of a Two-Dimensional Array with a Parallel Feed for 
N, = N, = 64, 6, = 0°, 9, = 0°. Contour Plot at 20 dB Level 
(Approximate Method). 
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In Band RCS of a two-dimensional array - Contour plot 
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Figure 15. Inband RCS of a Two-Dimensional Array with a Parallel Feed for 
N, = N, = 64, 8, = 45°, @, = 0°. Contour Plot at 20 dB Level 
(Approximate Method). 
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In Band RCS of a two-dimensional array - Contour plot 
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Figure 16. Inband RCS of a Two-Dimensional Array with a Parallel Feed for 
N, = Ny = 64, 8, = 45°, o, = 45°. Contour Plot at 20 dB Level 
(Approximate Method). 
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B. RCS DATA FOR THE RIGOROUS SOLUTION 

For the rigorous method, pattern plots of the RCS of linear arrays are 
obtained by using the FORTRAN program developed in [Ref. 6]. This program 
computes the RCS per squared wavelength in dB for any monostatic angle 0. 
Array parameters include the number of radiating elements V (N has to be a 
multiple of 8). scan angle 8, , and electrical path length between devices wy, (in 
radians). The output is written to MATLAB files that can be used to plot the RCS 
patterns. 

Figure 17 shows the RCS values of a linear array for N = 64 with no 
scanning (8, = 0) and path lengths of w, = 0 radians. As in the approximate 
method, the specular lobe at @ = 0, coupler lobes, and Bragg lobes at 8 = +85 
degrees are evident. The same array is examined in Figures 18 and 19 with wy, = 
n/4 and 1/2, respectively. A comparison of Figures 17 through 19 shows that the 
lobes have the same position but differ in amplitude. Specifically. for yw, = 7/4 the 
specular and coupler lobes have bigger amplitude than for y, = 0. However, for 
Wo = 1/4, the Bragg lobes have smaller amplitude than for y, = 0. The variation 
is due to the beating of mismatches in the feed. The addition or cancellation 
depends on the line lengths connecting the devices, wy. This effect becomes more 
complicated as the beam is scanned (Figure 20). Some lobes may disappear at 


some scan angles because the phase shifters introduce a shift that cause complete 
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Figure 17. Inband RCS of a Linear Array with a Parallel Feed for N = 64, 
8, = 0°, vw. = 0 (Rigorous Method). 
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Figure 18. Inband RCS of a Linear Array with a Parallel Feed for N = 64. 
8, = 0°. wy = 2/4 (Rigorous Method). 
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Figure 19. Inband RCS of a Linear Array with a Parallel Feed for N = 64, 
. = O°. wy = 2/2 (Rigorous Method). 
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Figure 20. Inband RCS of a Linear Array with a Parallel Feed for N = 64, 
6, = 45°, w, = 7/4 (Rigorous Method). 
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cancellation of two scattering contributions. This is not predicted by the 


approximate method because it sums the individual contributions noncoherently. 


C. COMPARISON SUMMARY 

The rigorous and approximate results have been presented for broadside and 
scanned linear arrays. In both cases the specular lobes have almost the same 
magnitude (within about 1dB) but the coupler lobes vary about 3 dB. This is 
attributed to the noncoherent summation in the approximation. When the beam is 
scanned, both methods predict the proper lobe locations. unless complete 
cancellation occurs. 

From a practical point of view, the approximate results are very close to the 
rigorous. One major difference is the computation times. For V = 64, the rigorous 
method takes about eight times longer than the approximate method and increases 
dramatically when N increases (for V = 128 over 12 times). A more important 
difference is that if the number of levels of couplers is increased. the scattering 
equations must be completely rewritten and programmed. For the approximate 


method, more terms only need to be added to equation (33). 


V. CONCLUSIONS 


An approximate scattering model for arrays with parallel feed networks has 
been presented. Calculations for several cases were compared to a rigorous method 
which includes all the interactions between the feed devices and aperture. The 
approximate method was in good agreement with the rigorous method in predicting 
RCS lobe positions, heights, and behavior with scanning. 

There are several advantages to the approximate approach. First, it is 
computationally efficient. allowing two-dimensional contours to be generated in 
minutes. Second, it can be easily extended to an arbitrary number of elements and 
coupler levels. The disadvantage is that a noncoherent addition of terms does not 
predict total cancellation conditions. However, most RCS designers are primarily 
concerned with the "worst case" conditions for highest RCS, and in this sense the 
approximate method is sufficient. 

Future efforts should be directed at increasing the number of coupler levels. 
and adding a coupling network in the y dimension of the two-dimensional array. 


Also. methods of reducing the inband RCS should be investigated. 
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APPENDIX A 
MATLAB PROGRAM FOR LINEAR ARRAYS 


% Phased arrays 
% Written by V.FLOKAS - 22 APR.1994 
% Inband RCS of a linear array with parallel feed 
clear 
clg 
thetas=input (’ Enter scanned angle in degrees thetas= ?) 
% Scanning at Phis=0 
this=thetas*pi/180; 
c=3e8; 
f=input(’Enter operating frequency in Hz, f= ’) 
lamda=c/f; 
k=2*pi/lamda; 
4=0.5*lamda; 

-chitk*d*sin(this) ; 
theta=linspace(-89,89, 660) ; 
thi=theta*pi/180; 

-alpha=k*d*sin(thi) ; 
jeta=alph4-chi*ones(1,length(alpha)); 
1=0.5*lamda; 

Ae=d*1; 

N=input(’Enter the number of elements N= ’) 
A=N*Ae; 
r=0.2; 
tesqrt(1-(r°2)); 
ai=(sin(N*alpha) ) ./(N*sin(alpha)) ; 
hi=find(isnan(a1)); 
ai(Chi)=ones(size(hi)); 
a2=(sin(N*jeta))./(N*sin(jeta)); 
h2=find(isnan(a2)); 
a2(h2)=ones(size(h2)); 
a3=(sin(N*jeta)) ./((N/2)*sin(2*jeta)); 
h3=find(isnan(a3)); 
a3(h3)=ones(size(h3)); 
a4=(sin(N*jeta)) ./((N/4)*sin(4*jeta)); 
h4=find(isnan(a4)); 
a4(h4)=ones(size(h4)); 
sigma=(4*pi*(A/lamda)~2)*((cos(thi)).72).*((r72)*(a1.72)+( (172) *(t74) * 

(ad. °2))+((t7°8)* (272) *(a2.72))+((67°8) * (272) * ((cos(jeta/2)).°4).* 
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(a3.72))+((t78) *(r72)*((sin(jeta/2)) .74) .*(a3.°2))+ 
((t78)* (272) *((cos(jeta/2)) .~4) .*((cos(jeta)).74) .*(a4.72))+ 
((t78)*(r72)*((cos(jeta/2)).°4).*((sin(jeta)).*4) .*(a4.°2))); 
sig=10*logi0(sigma/(lamda“~2)) ; 
figure(1) 
plot (theta,sig),grid 
xlabel(’Monostatic Angle (deg)’) 
ylabel(’sigma/lamda*2 (dB)’) 
title([’In Band RCS of a linear array with parallel feed N= ’ ,num2str(N),’ 
and thetas=’ ,num2str(thetas), J) 
gtext(’Approximate Method’) 
axis([-90 90 -50 40]) 
pause 
sec=(4*pi*(A/lamda) ~2)*((cos(thi)) .72) .*(((¢78) *(r72)*((cos(jeta/2)) .°4) .* 
((cos(jeta)) .74).*(a4.°72))+( (478) *(r°2) *((cos(jeta/2)) .74) .*((sin(jeta)) .74). 
*(a4.°72))); 
s=10*log10(sec/(lamda*2)) ; 
figure(2) 
plot(theta,s),grid,xlabel(’Monostatic Angle (deg)’),ylabel(’sigma/lamda”2 
(dB) ?) 
title([’In Band RCS of a linear array with Parallel feed N= ’,num2str(N),’ 
and thetas=’ ,num2str(thetas), ]) 
gtext(’Second level of couplers’) 
axis([-90 90 -50 40]) 
pause 
fir=(4*pi*(A/lamda) ~2)*((cos(thi)) .*2).*((t78)*(r72)*((cos(jeta/2)).74) .* 
((cos(jeta)) .74) .*(a4.72))+ 
((t78) * (172) *((cos(jeta/2)) .74) .*((sin(jeta) .74).*(a4.°2))); 
si=10*log10(fir/(lamda*2)) : 
figure (3) 
plot(theta,s1),grid,xlabel(’Monostatic Angle (deg)’),ylabel(’sigma/lamda*2 
(dB) ?) 
title([’In Band RCS of a linear array with parallel feed N= ’ ,num2str(N),’ 
and thetas=’ ,num2str(thetas), ]) 
gtext(’First level of couplers’) 
axis([-90 90 -50 40]) 


APPENDIX B 
MATLAB PROGRAM FOR PLANAR ARRAYS 


% Phased arrays 

% Written by V.FLOKAS - 22 APR.1994 

% Part (a) - Linear plot for 2-d array 

Clear 

clg 

thetas=input(’Enter scanned angle in degrees thetas= ’) 
this=thetas*pi/180; 

Phis=input(’Enter scanned angle in degrees Phis= ’) 
phis=Phis*pi/180; 

c=3e8; 

f=input(’Enter operating frequency in Hz= ’) 
lamda=c/f; 

k=2*pi/lamda; 

dx=0 .5*lamda; 

dy=dx; 

chiox=k*dx*sin (this) *cos(phis) ; 
chioy=k*dy*sin(this)*sin(phis) ; 
theta=linspace(-89,89,460) ; 

thi=theta*pi/180; 

PHI=Phis; % IN MONOSTATIC RCS (PHI CUT = PHI SCANNED) 
phi=PH1*pi/180; 
u=sin(thi).*(ones(1,length(thi))*cos(phi)); 
vesin(thi).*(ones(1,length(thi))+*sin(phi)); 
alpha=k*dx*u; 

beta=k*dy*v ; 
jetax=alpha-chiox*ones(1,length(alpha)) ; 
jetay=beta-chioy*ones(1,length(beta)); 
1=0.5*lamda; 

Nx=input (’Enter the number of elements Nx= ’) 
Ny=Nx; 

A=Nx*Ny*dx*dy ; 

r=0.2; 

t=sqrt(1-(r72)); 
ai=(sin(Nx*alpha))./(Nx*sin(alpha)) ; 
hi=find(isnan(al1)); 

ai(hi)=ones(size(h1)); 

a2=(sin(Nx*jetax)) ./(Nx*sin(jetax)); 
h2=find(isnan(a2)); 
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a2(h2)=ones(size(h2)); 
a3=(sin(Nx*jetax)) ./((Nx/2)*sin(2*jetax)); 
h3=find(isnan(a3)); 
a3(h3)=ones(size(h3)); 
a4=(sin(Ny*beta)) ./(Ny*sin(beta)); 
h4=find(isnan(a4)); 
a4(h4)=ones (size(h4)); 
a5=(sin(Ny*jetay))./(Ny*sin(jetay)) ; 
h5=find(isnan(a5)); 
a5(h5)=ones (size(h5)); 
a6=(sin(Nx*jetax))./((Nx/4)*sin(4*jetax)); 
hé=find(isnan(a6)); 
a6(h6)=ones(size(h6)); 
sigma=(4*pi*(A/lamda)~2)*((cos(thi)) .°2).*((r72)*((al.*a4@) .72)+( (172) * 
(t74)* CCad.*ad) .72))+0(478) * (1°72) *( (a2. *a5) 72) )+( (+78) * (172) * 
((cos(jetax/2)).74) .* 
((a3.*a5) .72))+((+78) *(r72) *((sin(jetax/2)).74) .*((a3.*a5) .72))+ 
((t 78) *(r°2)*((cos(jetax/2)) .74) .*((cos(jetax)) .74) .*((a6.*aS).72))+ 
( (678) *(r°2) *((cos(jetax/2)) .-4) .*((sin(jetax)).74) .*((a6.#a5).72))); 
sig=10*1logi0(sigma/(lamda*2)); 
% Contribution from first level of couplers (sum and difference arms) 
sigmal=(4*pi*(A/lamda)~2)*((cos(thi)).°2) .*(((t78)*(r72)* 
((cos(jetax/2)).*4). 
*( (a3. *a5) .72))+((t78) * (172) *C(sin(jetax/2)) .74).*((a3.*a5).72))); 
sigi=10*logi0(sigmai/(lamda”2)); 
% Contribution from second level of couplers 
sigma2=(4*pi*(A/lamda)~2)*((cos(thi)) .72).*(((t78) *(r°2)*((cos(jetax/2)). 
“4).*((cos(jetax)).°4).*((a6.*a5) .72))+( (+78) *(r°2)*((cos(jetax/2)).74). 
*((sin(jetax)) .74).*((a6.*a5).°72))); 
sig2=10*logi0(sigma2/(lamda*2)) ; 
figure(1) 
plot (theta,sig) ,grid 
axis([-90 90 -30 80]) 
xlabel(’Monostatic Angle (deg)’) 
ylabel (’sigma/lamda*2 (dB)’) 
title(’In Band RCS of a planar 2-D array with parallel feed’) 
gtext([’Approximate Method ~- Nx=’,num2str(Nx),’=Ny=N and thetas=’, 
num2str(thetas) 
gtext([’Phis=’ ,num2str(Phis), ]) 
pause 
figure (2) 


wr 


plot (theta,sigi),grid 
axis([-90 90 ~30 80]) 
xlabel(’Monostatic Angle (deg)’) 
ylabel(’sigma/lamda*2 (dB)’) 
title(’In Band RCS of a planar 2-D array with parallel feed - 
First Level of Couplers’) 
gtext([’Approximate Method - Nx=’ ,num2str(Nx),’=Ny=N and thetas=’ 
shum2str(thetas), ]) 
gtext((’Phis=’ ,num2str(Phis), ]) 
pause 
figure(3) 
plot (theta, sig2) ,grid 
axis([-90 90 -30 80]) 
xlabel(’Monostatic Angle (deg)’) 
ylabel(’sigma/lamda“2 (dB)’) 
title(’In Band RCS of a planar 2-D array with parallel feed - 
Second Level of Couplers’) 
gtext([’ Approximate Method - Nx=’ ,num2str(Nx),’=Ny=N and thetas= 
> num2str(thetas), ]) 
gtext([’Phis=’ ,num2str(Phis), ]) 
pause 
% Part(b) - Contour plot for the same array 
clear thi phiuv 
du=0.005; 
dv=0 .005; 
nu=fix(2/du)+1; 
nv=fix(2/dv)+1; 
for i=i:nu 
u(i)#-1+(i-i)*du; 
for j=i:nv 
v(j)8-1+(j-1) ¥dv; 
ec=(u(i))72+(v(j))°2; 
eth(i,j)=1- (Cui) °2) +(v(j)72)); 
if (ce>1) 
cth(i,j)=0; 
end 
end 
end 
[x,y]=meshgrid(u,v); 
alpha=k*dx*x; 
beta=k*dy*y ; 
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jetax=alpha-chiox*ones(length(alpha) , length(alpha)); 
jetay=beta-chioy*ones(length(beta) , length(beta)) ; 
ai=(sin(Nx*alpha)) ./(Nx*sin(alpha)) ; 

zi=find(isnan(al)); 

ai(zi)=ones(size(z1)); 

a2=(sin(Nx*jetax))./(Nx*sin(jetax)); 

z2=find(isnan(a2)); 

a2(z2)=ones(size(z2)); 
a3=(sin(Nx*jetax))./((Nx/2)*sin(2*jetax)); 

23=find(isnan(a3)); 

a3(z3)=ones(size(z3)); 

a4=(sin(Ny*beta)) ./(Ny*sin(beta)); 

z4=find(isnan(a4)); 

a4(z4)=ones(size(z4)); 

a5=(sin(Ny*jetay))./(Ny*sin(jetay)) ; 

25=find(isnan(a5)); 

a5(z5)=ones(size(z5)); 
a6=(sin(Nx*jetax))./((Nx/4)*sin(4*jetax)); 

z6=find(isnan(a6)); 

a6(z6)=ones(size(z6)); 

q=(4*pi*(A/lamda) °2)*(cth) ; 
n=(r72)*((at.#a4) .72)+(r72) *(t74)*( (al .*a4d) .72)+(t78)* (172) * 
((a2.*a5) .72)4+(t78) * (272) *( (cos(jetax/2)).7"4).* 

((a3.*a5) .72)+(¢78)*(r72)*((sin(jetax/2)) .-4).*((a3.*a5).72)+ 
(t°8)*(r72) *((cos(jetax/2)) .74).*((cos(jetax)) .-4) .*((a6.*a5) .72)+ 
(478) * (172) *((cos(jetax/2)).74) .*((sin(jetax)) .74) .*((a6.*a5) .°2); 

sigma=q.*n; 

Sigma=abs (sigma) ; 
sig=10*logi0((sigma/(lamda“2))+eps*ones(length(sigma) , length(sigma))); 
lev=[80 20]; 

figure (4) 

axis(’ square’) 

contour(sig,lev,u,v),grid 

axis(’square’) 

xlabel(?u’) ,ylabel(’v’) 

title(’In Band RCS of a two-dimensional array - Contour plot’) 
gtext([’ Approximate Method Nx=’ ,num2str(Nx) ,’=Ny=N and thetas=’ 
pnum2str(thetas) ,] 

gtext([’?Phis=’ ,num2str(Phis), ]) 


to 
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